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We report on the high statistics two-pion correlation functions from pp collisions at = 0.9 TeV and 
i/s = 7 TeV, measured by the ALICE experiment at the Large Hadron Collider. The correlation functions as 
well as the extracted source radii scale with event multiplicity and pair momentum. When analyzed in the same 
multiplicity and pair transverse momentum range, the correlation is similar at the two collision energies. A three- 
dimensional femtoscopic analysis shows an increase of the emission zone with increasing event multiplicity as 
well as decreasing homogeneity lengths with increasing transverse momentum. The latter trend gets more 
pronounced as multiplicity increases. This suggests the development of space-momentum correlations, at least 
for collisions producing a high multiplicity of particles. We consider these trends in the context of previous 
femtoscopic studies in high-energy hadron and heavy-ion collisions, and discuss possible underlying physics 
mechanisms. Detailed analysis of the correlation reveals an exponential shape in the outward and longitudinal 
directions, while the sideward remains a Gaussian. This is interpreted as a result of a significant contribution 
of strongly decaying resonances to the emission region shape. Significant non-femtoscopic correlations are 
observed, and are argued to be the consequence of "mini-jet"-like structures extending to low pj. They are well 
reproduced by the Monte-Carlo generators and seen also in 71+71^ correlations. 

PACS numbers: 25.75.-q, 25.75.Gz, 25.70.Pq 
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Proton-proton collisions at y/s = 0.9 TeV and y/s — 7 TeV 
have been recorded by A Large Ion Collider Experiment (AL- 
ICE) at the Large Hadron Collider (LHC) at CERN in 2010. 
These collisions provide a unique opportunity to probe Quan- 
tum Chromodynamics (QCD) in the new energy regime. The 
distinguishing feature of QCD is the mechanism of color con- 
finement, the physics of which is not fully understood, due, 
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in part, to its theoretical intractability iQl]. The confinement 
mechanism has a physical scale of the order of the proton ra- 
dius and is especially important at low momentum. The study 
presented in this work aims to measure the space-time extent 
of the source on this scale. 

Two-pion correlations at low relative momentum were first 
shown to be sensitive to the spatial scale of the emitting source 
in p + p collisions by G. Goldhaber, S. Goldhaber, W. Lee 
and A. Pais SOyears ago |2j. Since then, they were stud- 
ied in e+ + yfl, hadron- and lepton-hadron iQ], and heavy 
ion collisions. Especially in the latter case, two-particle 
femtoscopy has been developed into a precision tool to probe 
the dynamically-generated geometry structure of the emitting 
system. In particular, a sharp phase transition between the 
color-deconfined and confined states was precluded by the ob- 
servation of short timescales, and femtoscopic measurement 
of bulk collective flow proved that a strongly self-interacting 
system was created in the collision i^l?!]- 

Femtoscopy in heavy-ion collisions is believed to be under- 
stood in some detail, see e.g. jsj]. The spatial scales grow nat- 
urally with the multiplicity of the event. Strong hydrodynami- 
cal collective flow in the longitudinal and transverse directions 
is revealed by dynamical dependencies of femtoscopic scales. 
The main puzzling aspect of the data is the relative energy 
independence of the results of the measurements. 

To some extent, Bose-Einstein correlations in particle 
physics were initially of interest only as a source of sys- 
tematic uncertainty in the determination of the W boson 
mass (sl. But overviews (3, 4, !'9!] of femtoscopic measure- 
ments in hadron- and lepton-induced collisions reveal system- 
atics surprisingly similar to those mentioned above for heavy- 
ion collisions. Moreover, in the first direct comparison of fem- 
toscopy in heavy-ion collisions at Relativistic Heavy-Ion Col- 
lider (RHIC) and proton collisions in the same apparatus an 
essentially identical multiplicity- and momentum-dependence 
was reported in the two systems ifioll . However, the multiplic- 
ities at which the femtoscopic measurement in pp collisions at 
RHIC was made were still significantly smaller than those in 
even the most peripheral heavy-ion collisions. In this work we 
are, for the first time, able to compare femtoscopic radii mea- 
sured in pp and heavy-ion collisions at comparable event mul- 
tiplicities. At these multiplicities the observed correlations 
may be influenced by jets jl 1 I J while other studies suggest that 
a system behaving collectively may be created lfT2ll . 

In our previous work ifisll we reported that a multiplicity 
integrated measurement does not show any pair momentum 
dependence of the /?,,„, radius measured in the Pair Rest Frame 
(PRE). Similar analysis from the CMS collaboration ifTill also 
mentions that no momentum dependence was observed. How- 
ever the analysis in two multiplicity ranges suggested that mo- 
mentum dependence may change with multiplicity, although 
any strong conclusions were precluded by limited statistics. In 
this work we explore this dependence by using high statistics 
data and more multiplicity ranges. It enabled us to perform the 
three-dimensional analysis in the Longitudinally Co-Moving 
System (LCMS), where the pair momentum along the beam 
vanishes. 

The paper is organized as follows: in Section|II]we describe 



the ALICE experimental setup and data taking conditions for 
the sample used in this work. In Section|lIl]we present the cor- 
relation measurement and characterize the correlation func- 
tions themselves. In Section lTV Al we show the main results of 
this work: the three-dimensional radii extracted from the data. 
We discuss various observed features and compare the results 
to other experiments. In Section |V] we show, for complete- 
ness, the one-dimensional 7?,,,,, analysis. Finally in SectionlVll 
we summarize our results. All the numerical values can be 
obtained from the Durham Reaction Database ifisll . 



II. ALICE DATA TAKING 

In this study we report on the analysis of pp collisions 
recorded by the ALICE experiment during the 2010 run of the 
LHC. Approximately 8 million events, triggered by a mini- 
mum bias trigger at the injection energy of y/s = 0.9 TeV, and 
100 million events with similar trigger at the maximum LHC 
energy to date, y/s = 7 TeV, were analyzed in this work. 

The ALICE Time Projection Chamber (TPC) iH] was used 
to record charged particle tracks as they left ionization trails 
in the Ne — CO2 gas. The ionization drifts up to 2.5 m 
from the central electrode to the end-caps to be measured on 
159 padrows, which are grouped into 18 sectors; the posi- 
tion at which the track crossed the padrow was determined 
with resolutions of 2 mm and 3 mm in the drift and transverse 
directions, respectively. The momentum resolution is ~ 1 % 
for pions with pj = 0.5 GeV/c. The ALICE Inner Tracking 
System (ITS) was also used for tracking. It consists of six sili- 
con layers, two innermost Silicon Pixel Detector (SPD) layers, 
two Silicon Drift Detector (SDD) layers, and two outer Sili- 
con Strip Detector (SSD) layers, which provide up to six space 
points for each track. The tracks used in this analysis were 
reconstructed using the information from both the TPC and 
the ITS, such tracks were also used to reconstruct the primary 
vertex of the collision. For details of this procedure and its ef- 
ficiency see UtIi . The forward scintillator detectors VZERO 
are placed along the beam line at +3 m and —0.9 m from the 
nominal interaction point. They cover a region 2.8 < r| < 5.1 
and —3.7 < r| < —1.7 respectively. They were used in the 
minimum bias trigger and their timing signal was used to re- 
ject the beam-gas and beam-halo collisions. 

The minimum bias trigger required a signal in either of the 
two VZERO counters or one of the two inner layers of the Sil- 
icon Pixel Detector (SPD). Within this sample, we selected 
events based on the measured charged-particle multiplicity 
within the pseudorapidity range |r|| < 1.2. Events were re- 
quired to have a primary vertex within 1 mm of the beam line, 
and 10 cm of the center of the 5 m-long TPC. This provides 
almost uniform acceptance for particles with |r| | < 1 for all 
events in the sample. It decreases for 1.0 < |r|| < 1.2. In ad- 
dition, we require events to have at least one charged particle 
reconstructed within |r|| < 1.2. 

The minimum number of clusters associated to the track 
in the TPC is 70 (out of the maximum of 159) and 2 in the 
ITS (out of the maximum of 6). The quality of the track is 
determined by the X^/N value for the Kalman fit to the re- 
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constructed position of the TPC clusters (A^ is the number of 
clusters attached to the track); the track is rejected if the value 
is larger than 4.0 (2 degrees of freedom per cluster). Tracks 
with |r|| < 1.2 are taken for the analysis. The pj of accepted 
particles has a lower limit of 0.13 GeV/c, because tracks with 
lower pj do not cross enough padrows in the TPC. The effi- 
ciency of particle reconstruction is about 50% at this lowest 
limit and then quickly increases and reaches a stable value of 
approximately 80% for pj>Q.2 Ge V/c. In order to reduce the 
number of secondary particles in our sample, we require the 
track to project back to the primary interaction vertex within 
0.018 + 0.035/?j cm in the transverse plane and 0.3 cm in 
the longitudinal direction (so-called Distance of Closest Ap- 
proach or DCA selection). 

ALICE provides an excellent particle identification capabil- 
ity, through the combination of the measurement of the spe- 
cific ionization {dE /dx) in the TPC and the ITS and the timing 
signals in the ALICE Time Of Flight (TOE). In the momentum 
range covered here (0.13 Ge V/c to 0.7 Ge V/c) pions constitute 
the majority of particles. We use only the TPC measurement 
for Particle IDentification (PID) in this work, as the other de- 
tectors offer significant improvement at higher pj than used 
here. This PID procedure results in a small contamination of 
the pion sample by electrons at pj < 0.2 GeV/c and kaons at 
Pj > 0.65 GeV/c. Allowing other particles into our sample 
has only a minor effect of lowering the strength of the correla- 
tion (the A, parameter), while it does not affect the femtoscopic 
radius, so we do not correct for it explicitly. The amount of 
electron contamination is less than 5%. 



III. CORRELATION FUNCTION MEASUREMENT 

Experimentally, the two-particle correlation function is de- 
fined as the ratio C(q) =A(q)/Z?(q), where A(q) is the 
measured two-pion distribution of pair momentum difference 
q = P2 — pi, and B (q) is a similar distribution formed by us- 
ing pairs of particles from different events ifTsIl . 

The size of the data sample used for this analysis al- 
lowed for a highly differential measurement. In order to ad- 
dress the physics topics mentioned in the introduction, the 
analysis was performed simultaneously as a function of the 
total event multiplicity Nd, and pair transverse momentum 
kj = \pt.\ + Pt.iI /2. Eor the multiplicity determination we 
counted the tracks reconstructed simultaneously in the ITS 
and the TPC, plus the tracks reconstructed only in the ITS 
in case the track was outside of the TPC r| acceptance. The 
total number of events accepted after applying the selection 
criteria in the = 7 TeV sample was 60 x 10^ and in the 
y/s = 0.9 TeV sample it was 4.42 x lO*". We divided the full 
multiplicity range into eight and four ranges for the two ener- 
gies respectively in such a way that the like-charge pion pair 
multiplicity in each of them was comparable. TableUgives (a) 
values for the range of raw charged particle multiplicity that 
was used to categorize the event, (b) the corresponding mean 
charge particle density {dNch/dr[) as well as (c) number of 
events and (d) the number of identical pion pairs in each range. 
The femtoscopic measurement requires the events to have at 



TABLE I. Multiplicity selection for the analyzed sample. Uncor- 
rected Nch in |r|| < 1.2, (d/Vf./,/dr|) |A'd,>l (see text for the definition), 
number of events and number of identical pion pairs in each range 
are given. 



Bin 


Nch 


{dN,,,/dn) 


N^i,>l No. events x 10^ 


No. pairs x 10^ 


V^ = 0.9TeV 


1 


1-11 


2.7 


3.1 


8.8 


2 


12-16 


7.0 


0.685 


8.6 


3 


17-22 


9.7 


0.388 


9.5 


4 


23-80 


14.6 


0.237 


12.9 


^ = 7 TeV 


1 


1-11 


3.2 


314 


48.7 


2 


12-16 


7.4 


9.2 


65.0 


3 


17-22 


10.4 


7.4 


105.7 


4 


23-29 


13.6 


4.8 


120.5 


5 


30-36 


17.1 


3.0 


116.3 


6 


37^4 


20.2 


2.0 


115.6 


7 


45-57 


24.2 


1.3 


114.5 


8 


58-149 


31.1 


0.72 


108.8 
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'^outside' 



O 
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out.side.long 



FIG. 1. Projections of the 3D Cartesian representations of the corre- 
lation functions onto the qou,, qside^ and qiong axes for pairs with 
0.2 < kj < 0.3 GeV/c, for three multiplicity ranges. To project 
onto one ^-component, the others are integrated over the range 
0-0.16 GeV/c. 



' In fact the coirelation signal is constructed from events having at least two 
same-charge pions (a pair). The one-pion events do contribute to the mixed 
background. 
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sample. We denote this value as (dA^t./,/dr|) |a?^,,>i; its typical 
uncertainty is 10%. We note that for the the lowest multi- 
plicity this charged particle density is biased towards higher 
values with respect to the full sample of inelastic events. 

The pair momentum ranges used in the analysis 
were (0.13,0.2), (0.2,0.3), (0.3,0.4), (0.4,0.5), (0.5,0.6), 
(0.6,0.7) GeV/c. 



A. Correlation function representations 

The correlations are measured as a function of pair relative 
momentum four-vector q. We deal with pions, so the masses 
of the particles are fixed - in this case q reduces to a vec- 
tor: q. The one-dimensional analysis is performed versus the 
magnitude of the invariant momentum difference ^,„v = \q\, 
in PRF. The large available statistics for this work allowed 
us to perform a detailed analysis also for the 3D functions. 
In forming them, we calculate the momentum difference in 
LCMS and decompose this qtcMS according to the Bertsch- 
Pratt lfl9[ I20I1 "out-side-long" (sometimes indicated by o, s, 
and I subscripts) parametrization. Here, qiong is parallel to 
the beam, qow is parallel to the pair transverse momentum, 
and qside is perpendicular to qtong and qom- If one wishes to 
compare the radii measured in LCMS to /?,„,, one needs to 
multiply one of the transverse radii in LCMS (the one along 
the pair transverse momentum) by the Lorentz y correspond- 
ing to the pair transverse velocity, and then average the three 
radii. Therefore an /?,,„, constant with momentum is consistent 
with the radii in LCMS decreasing with momentum. Figure[T| 
shows one-dimensional projections of the 3-dimensional cor- 
relation function C{qout,qside,qiong) onto the ^„„,, q^ide, and 
qiong axes, for 7t+ pairs from one of the multiplicity/fcx ranges 
from the ^/s = 7 Te V sample. The function is normalized with 
a factor that is a result of the fit (the details of the procedure 
are described in Sec. lIIIDl ): unity means no correlation. 

The 1 -dimensional projections, shown in Fig. [T] present 
a limited view of the 3-dimensional structure of the correla- 
tion function. It is increasingly common to represent correla- 
tion functions in a harmonic analysis ll2Ti - l23ll : this provides a 
more complete representation of the 3-dimensional structure 
of the correlation, a better diagnostic of non-femtoscopic cor- 
relations 12211. and a more direct relation to the shape of the 
source 112411 . The moments of the Spherical Harmonic (SH) 
decomposition are given by 



1 



^/^^/(cose)c(|^|,e,^)y/"(0,(^). (i) 



Here, the out-side-long space is mapped onto Euler angles in 
which qiong ~ |^|cos9 and ^„„, = |^| sin9cos(|). For pairs of 
identical particles in collider experiments done with symmet- 
rical beams, including the analysis in this work, the odd / and 
the imaginary and odd m components for even / vanish. The 
first three non-vanishing moments, which capture essentially 
all of the 3-dimensional structure, are then Cq, Cj, and C^- 
These are shown in Fig. |2l The components for / > 4 rep- 
resent the fine details of the correlation structure and are not 




■"LCMS 



1.0 

(GeV/c) 



FIG. 2. Moments of the SH decomposition of the correlation func- 
tions for pairs with 0.2 < < 0.3 GeV/c, for three multiplicity 
ranges. 



analyzed in this work. 

The Cq is the angle-averaged component. It captures the 
general shape of the correlation. The width of the peak near 
(7 = is inversely proportional to the overall femtoscopic size 
of the system. The component is the correlation weighed 
with the cos^(9). If it differs from 0, it signifies that the longi- 
tudinal and transverse sizes of the emission region differ. The 
C2 is weighed with cos^ ((])). If it differs from 0, it signals that 
the outward and sideward sizes differ. The correlation func- 
tion is normalized to the number of pairs in the background 
divided by the number of pairs in the signal. 



B. Measured correlations 

In Figs. [U and |2] we show selected correlations to illus- 
trate how they depend on multiplicity. This is done for kj 
of (0.2,0.3) GeV/c; the behavior in other kj ranges and at 
the lower collision energy is qualitatively the same. The nar- 
rowing of the correlation peak with increasing multiplicity is 
apparent, corresponding to the increase of the size of the emit- 
ting region. The behavior of the correlation function at large q 
is also changing, the low multiplicity baseline is not flat, goes 
below 1 .0 around q= \ GeV/c and then rises again at larger 
q, for higher multiplicities the background becomes flatter at 
large q. In Cartesian representation shown in Fig.[Tl areas with 
no data points (acceptance holes) are seen in qout projection 
near q = 0.5 GeV/c and in qiong above 0.6 GeV/c. Since qiong 
is proportional to the difference of longitudinal momenta, its 
value is limited due to r| acceptance. In the out direction the 
hole appears due to a combination of lower pi cut-off and the 
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FIG. 3. Moments of the SH decomposition of the correlation func- 
tions for events with 17 < Nch < 22, for three kj ranges. 



selected kr range. It can be simply understood as follows: For 
the projection in the upper panel of Fig.[Tl we take the value 
of 

qside ^nd qiong small. The value of qside is proportional 
to the azimuthal angle difference, while qioni; is proportional 
to polar angle difference. For qside, qumg = 0, qout is simply 
Pt.2 — Pt.i and kj is {pjA + Pt.i) /2, where pj is no longer a 
2-vector, but just a scalar. The particles are either fully aligned 
(both pj's are positive or both are negative) or back-to-back 
(one pt is positive, the other negative). When we combine 
the lower pj cut-off \pt\ > 0.13 GeV/c and the kj selection 
0.2 < kj < 0.3, it can be shown that some range of the qout 
values is excluded. This range will depend on the selec- 
tion. 

The kj dependence of the correlation function is shown in 
Figs.[3]and|4] for multiplicity 17 < Nd, < 22. The behavior in 
other multiplicity ranges and at lower energy is qualitatively 
similar (except the lowest multiplicity bin where the behav- 
ior is more complicated - see the discussion of the extracted 
radii in Sec. IIIIDj for details). We see a strong change of the 
correlation with kj, with two apparent effects. At low kj the 
correlation appears to be dominated by the femtoscopic effect 
at ^ < 0.3 GeV/c, and is flat at larger q. As kj grows, the fem- 
toscopic peak broadens (corresponding to a decrease in size of 
the emitting region). In addition, a wide structure, extending 
up to 1.0 GeV/c in q for the highest range, appears. We 
analyze this structure in further detail later in this work. We 
also see that, according to expectations, the acceptance holes 
in the out and long region move as we change the kj range. 

Figure|5]shows the example of the correlation function, for 
the same multiplicity /fcx range, for pp collisions at two col- 
lision energies. We note a similarity between the two func- 
tions; the same is seen for other A:t's and overlapping mul- 



FIG. 4. Projections of the 3D Cartesian representations of the corre- 
lation functions onto the qout, qside, ™d <ilong axes, for events with 
17 < A'c/, < 22, for three ranges. To project onto one ^-component, 
the others are integrated over the range — 0. 16 GeV/c. 
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tiplicity ranges. The similarity is not trivial: changing the 
multiplicity by 50%, as seen in Fig. |2]or kj by 30% as seen 
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in Fig. [3] has a stronger influence on the coiTelation function 
than changing the collision energy by an order of magnitude. 
We conclude that the main scaling variables for the correlation 
function are global event multiplicity and transverse momen- 
tum of the pair; the dependence on collision energy is small. 
The energy independence of the emission region size is the 
first important physics result of this work. We emphasize that 
it can be already drawn from the analysis of the correlation 
functions themselves, but we will also perform more quaUta- 
tive checks and discussions when we report the fitted emission 
region sizes in Section lTVl 



C. Non-femtoscopic correlation structures 

In Fig. [3] we noted the appearance of long-range structures 
in the correlation functions for large kj. If these were of fem- 
toscopic origin, they would correspond to an unusually small 
emission region size of 0.2 fm. We reported the observation of 
these structures in our previous analysis 11311 at ^/s = 0.9 TeV, 
where they were interpreted as non-femtoscopic correlations 
coming from "mini-jet" like structures at < 1 GeV/c. Here 
we further analyze this hypothesis. In Fig.|6]we show the com- 
parison of the correlation function at multiplicity 12 < Nd, < 
16 in an intermediate kj range, where the long-range corre- 
lations are apparent, to the Monte-Carlo (MC) calculation. 
The simulation used the PYTHIA generator 112511 . Perugia-0 
tune ||26|] as input and was propagated through the full sim- 
ulation of the ALICE detector 1 16-|. Then it was reconstructed 
and analyzed in exactly the same way as our real data, us- 



ing the same multiplicity and kj ranges. The MC calculation 
does not include the wave-function symmetrization for iden- 
tical particles; hence, the absence of the femtoscopic peak at 
low q is expected. In the angle-averaged Cq component a sig- 
nificant correlation structure is seen, up to 1 GeV/c, with a 
slope similar to the data outside of the peak at low q. Sim- 
ilarly, in the C? component a weak and wide correlation dip 
is seen around q = 0.5 GeV/c, which is also seen in the data. 
In MC, the correlation in C2 disappears at lower q, while for 
the data it extends to much lower q, exactly where the fem- 
toscopic peak is expected and seen in Cq. Our hypothesis is 
that both the long-range peak in and the dip in are of a 
"mini-jet" origin. They need to be taken into account when fit- 
ting the correlation function from data, so that the femtoscopic 
peak can be properly extracted and characterized. The calcu- 
lations was also carried out with a second Monte-Carlo, the 
PHOJET model lH^,!!!], and gave similar results. The differ- 
ences between the two models are reflected in the systematic 
error. 

In order to characterize the non-femtoscopic background 
we study in detail the correlation structure in the MC genera- 
tors, in exactly the same multiplicity /fcj ranges as used for data 
analysis. We see trends that are consistent with the "mini-jet" 
hypothesis. The correlation is small or non-existent for low 
Pj (first kj range) and it grows strongly with pj. In Fig.|2]we 
show this structure for selected multiplicity/A:T at both ener- 
gies. At the highest kj the effect has the magnitude of 0.3 at 
low q, comparable to the height of the femtoscopic peak. The 
appearance of these correlations is the main limiting factor in 
the analysis of the kj dependence. We tried to analyze the 
correlations at higher than 0.7 GeV/c but we were unable 
to obtain a meaningful femtoscopic result, because the "mini- 
jet" structure was dominating the coiTelation. The strength 
of the correlation decreases with growing multiplicity (as ex- 
pected), slower than 1 /M, so that it is still significant at the 
highest multiplicity. We studied other tunes of the PYTHIA 
model and found that the Perugia-0 tune reproduces the "mini- 
jet" correlation structures best, which is why it is our choice. 
Its limitation though is a relatively small multiplicity reach, 
smaller than the one observed in data. As a result the MC 
calculation for our highest multiplicity range is less reliable - 
this is reflected in the systematic error. 

Analyzing the shape of the underlying event correlation for 
identical particle pairs in MC is important; however, it does 
not ensure that the behavior of the correlation at very low q 
is reproduced well in MC. We compared the identical parti- 
cle MC and data in the large q region, where the femtoscopic 
effect is expected to disappear, and found them to be very 
similar in all multiplicity /fcx- However, if the "mini-jet" hy- 
pothesis is correct, the same phenomenon causes similar cor- 
relations to appear for non-identical pions. The magnitude is 
expected to be higher than for identical pions, because it is 
easier to produce an oppositely-charged pair from a fragment- 
ing "mini-jet" than it is to create an identically-charged pair, 
due to local charge conservation. Moreover, the femtoscopic 
effect for non-identical pions comes from the Coulomb inter- 
action only. It is limited to very low q, below 0.1 GeV/c. It is 
therefore possible to test the low-q behavior of the "mini-jet" 
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correlation with such correlations. In Fig. [8] we show the mea- 
sured 7t+7t^ correlation functions, in selected multiplicity/^x 
ranges, compared to the corresponding correlations from the 
same MC sample which was used to produce correlations in 
Fig. 13 The underlying event long-range correlation is well re- 
produced in the MC. We see some deviation in the lowest mul- 
tiplicity range, which is taken into account in the systematic 
error estimation. At larger multiplicities the strength of the 
correlation is well reproduced. By comparing the 3D function 
in SH we checked that the shape in 3D q space is also in agree- 
ment between data and MC. The magnitude for non-identical 
pions is slightly bigger than for identical pions, as expected. 
The femtoscopic Coulomb effect at < 0.1 GeV/c is also vis- 
ible. Another strong effect, even dominating at low multiplic- 
ity, are the peaks produced by the correlated pairs of pions 
coming from strong resonance decays. They do appear in the 
MC as well, but they are shifted and have different magnitude. 
This is the effect of the simplified treatment of resonance de- 
cays in PYTHIA, where phase space and final state rescattering 



are not taken into account. By analyzing some of the correla- 
tion functions in Fig [8] we were able to identify signals from 
at least the following decays: two-body p, /o, and /2 mesons 
decays, three-body co meson decay, and also possibly r| meson 
two-body decay. Some residual weak decay pairs, which 
are not removed by our DCA selection, can also be seen. All 
of these contribute through the full q range (0.0, 1.2) GeV/c. 
This fact, in addition to the stronger "mini-jet" contribution to 
non-identical (as compared to identical) correlations, makes 
the non-identical correlation not suitable for the background 
estimation for identical pion pairs. We also note that there 
appears to be very rich physics content in the analysis of res- 
onances decaying strongly in the 7t+7t^ channel; however, we 
leave the investigation of this topic for separate studies. 

The study of the 7t+7t^ correlations confirms that the MC 
generator of choice reproduces the underlying event structures 
also at low q. We found that they are adequately described by 
a Gaussian in LCMS for the Cq component. The dashed lines 
in Fig. |7] show the fit of this form to the correlation in MC. 
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The results of this fit, taken bin-by-bin for all multiplicity/A:T 
ranges, are the input to the fitting procedure described in Sec- 
tion llll DI Similarly, the observed correlation can be char- 
acterized well by a Gaussian, with the magnitude of —0.01 
or less and a peak around q ~ 0.5 GeV/c with a width of 
0.25 — 0.5 GeV/c. We proceed in the same way as for Cq; 
we fit the MC correlation structures with this functional form 
and take the results as fixed input parameters in the fitting of 
the measured correlations. 



D. Fitting tlie correlation function 

Having qualitatively analyzed the correlation functions 
themselves we move to the quantitative analysis. The femto- 
scopic part of the correlation function is defined theoretically 
via the Koonin-Pratt equation ll29ll30ll : 



(2) 



where q is the pair 3-momentum difference (the fourth com- 
ponent is not independent for pairs of identical pions since 
masses of particles are fixed), k is the pair total momentum, r 
is the pair space-time separation at the time when the second 
particle undergoes its last interaction, *P is the wave function 
of the pair, and S is the pair separation distribution. The aim in 
the quantitative analysis of the correlation function is to learn 
as much as possible about S from the analysis of the measured 
C. The correlation function C is, in the most general form, a 
6-dimensional object. We reduce the dimensionality to 3 by 
factorizing out the pair momentum k. We do not study the 
dependence on the longitudinal component of k in this work. 
The dependence on the transverse component of k is studied 
via the ki binning, introduced in Section|II] We assume that S 
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is independent of k inside each of the kj ranges. We also note 
that for identical pions the emission function 5 is a convolu- 
tion of two identical single particle emission functions . 

In order to perform the integral in Eq. Q we must postu- 
late the functional form of S or S\. We assume that S does 
not depend on q. The first analysis is performed with as 
a 3-dimensional ellipsoid with Gaussian density profile. This 
produces S which is also a Gaussian (with o larger by a factor 
of V2): 



(3) 



^'t'lmg P^*^'^ femtoscopic radii, also 



where Rg„ 

known as "HBT radii" or "homogeneity lengths", and 
and r/ are components of the pair separation vector. For iden- 
tical charged pions ^ should take into account the proper 
symmetrization, as well as Coulomb and strong interaction 
in the final state. In the case of the analysis shown in this 
work, with pions emitted from a region with the expected size 
not larger than 2 — 3 fm, the strong interaction contribution 
is relatively small and can be neglected ifsill . The influence 
of the Coulomb interaction is approximated with the Bowler- 
Sinyukov method. It assumes that the Coulomb part can be 
factorized out from ^ and integrated independently. There are 
well-known limitations to this approximation but they have 
minor influence for the analysis shown in this work. With 
these assumptions I' is a sum of two plane waves modified by 
a proper symmetrization. By putting Eq. (O into Eq. Q the 
integration can be done analytically and yields the quantum 
statistics-only correlation 



1+A.exp(-/?^ 



2 2 
1 lout ' 



pG 2 2 



?L.e)> (4) 



where A, is the fraction of correlated pairs for which both pi- 
ons were correctly identified. The 3-dimensional correlation 
function is then modified with the Bowler-Sinyukov formula 
to obtain the complete femtoscopic component of the correla- 
tion Cf: 



Cf{q)^{l-X)+XK{q„,) 



(5) 



l+exp( qlut-Rside Iside' 
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long 1U 
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where K is the Coulomb like-sign pion pair wave function 
squared averaged over the Gaussian source with a radius of 
1 fm. Changing this radius within the range of values mea- 
sured in this work has negligible effect on the extracted radii. 
Eq. (|5]l describes properly the femtoscopic part of the two 
pion correlation function. However, in the previous section 
we have shown that our experimental functions also contain 
other, non-femtoscopic correlations. We studied them in all 
multiplicity/^T ranges and found that they can be generally 
described by a combination of an angle-averaged Gaussian in 
LCMS plus a small Gaussian deviation in the component: 
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FIG. 9. Moments of the SH decomposition of the correlation func- 
tions for events with 23 < N^i, < 29 and pairs with 0.3 < < 
0.4 GeV/c. The dashed line shows the Gaussian fit, the dash-dotted 
line shows the background component, the dotted line show the fem- 
toscopic component. 
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(3cos2(e)-l), 



where A/,, A„, B/,, Bi„ and B„ are parameters. They are ob- 
tained, bin-by-bin, from the fit to the MC simulated correla- 
tion functions shown in Fig.|7] They are fixed in the procedure 
of fitting the data. The final functional form that is used for 
fitting is then: 



C{q 

out 1 Qside-, Qlong 

) = NCf{qout , qside , qiong)B {qout , qside , qiong), 

(7) 

where is the overall normalization. Projections of the 
Cartesian representation of the correlation functions, shown 
in Figs.[T]and|4l are normalized with this factor Function ^ 
is used to fit both the SH and Cartesian representation of the 
3D correlations. 

In Fig.|9]an example of the fit to one of our correlation func- 
tions is shown. The SH representation of the data is shown 
as points; the result of the fit is a black dashed line. The 
femtoscopic component is shown as a blue dotted line, the 
non-femtoscopic background as green dash-dotted line. The 
correlation function in this range has significant contribution 
from the background and is reasonably reproduced by the fit. 
At ^ < 0.1 GeV/c the fit misses the data points in Cq and C^; 
we will discuss what can be done to improve the agreement 
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in Section ITV CI In Fig. [TO] the same correlation is shown as 
projections of the 3D Cartesian representation. The other q 
components are integrated over the range of — 0.16 GeV/c. 
The fit, shown as lines, is similarly projected. In this plot 
the fit does not describe the shape of the correlation perfectly; 
however, the width is reasonably reproduced. 



rV. FIT RESULTS 
A. Results of the 3D Gaussian fits 

We fitted all 72 correlation functions (4h-8 multiplicity 
ranges for two energies times 6 ranges) with Eq. (|7|. We 
show the resulting femtoscopic radii in Fig.[TT]as a function of 
kj. The strength of the correlation X is relatively independent 
of kf, is 0.55 for the lowest multiplicity, decreases monoton- 
ically with multiplicity and reaches the value of 0.42 for the 
highest multiplicity range. The radii shown in the Fig.fTTIare 
the main results of this work. Let us now discuss many aspects 
of the data visible in this figure. 

Firstly, the comparison between the radii for two ener- 
gies, in the same multiplicity/fcx ranges reveals that they are 
universally similar, at all multiplicities, all ^t's and all di- 
rections. This confirms what we have already seen directly 
in the measured correlation functions. The comparison to 
^/s = 200 GeV pp collisions at RHIC is complicated by the 
fact that these data are not available in multiplicity ranges. 
The multiplicity reach at RHIC corresponds to a combination 
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of the first three multiplicity ranges in our study. No strong 
change is seen between the RHIC and LHC energies. It shows 
that the space-time characteristics of the soft particle produc- 
tion in pp collisions are only weakly dependent on collision 
energy in the range between 0.9 TeV to 7 TeV, if viewed in 
narrow multiplicity/fcx ranges. Obviously the ^/s = 7 TeV 
data have a higher multiplicity reach, so the minimum-bias 
(multiplicity/fcx integrated) correlation function for the two 
energies is different. 

Secondly, we analyze the slope of the ki dependence. Rf^^^^^ 
falls with kj at all multiplicities and both energies. and 
^^ide show an interesting behavior - at low multiplicity the ki 
dependence is flat for R^^^j^ and for it rises at low and 
then falls again. For higher multiplicities both transverse radii 
develop a negative slope as multiplicity increases. At high 
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FIG. 12. Gaussian radii vs event multiplicity, for ^/s = 0.9 TeV and 
7 TeV. Panel a) shows Rg,„, b) shows Rfi^^, c) shows Rfg^^g and d) 

shows Rout/Rfide ratio. Lines show linear fits to combined ^/s = 
0.9 TeV and i/s = 7 TeV points in each kj range. 



multiplicity the slope is bigger for while Rf^^^ grows uni- 
versally at all A:t's while developing a smaller negative slope. 
The difference in the evolution of shapes of R^^^^ and Rfj^^ is 
best seen in their ratio, shown in panel d) of Fig. (TT] At low 
multiplicities the ratio is close to 1.0, then it decreases mono- 
tonically with multiplicity. We note that a negative slope in 
Rg^if and R'^jj^ was universally observed in all heavy-ion mea- 
surements at RHIC energies and sometimes also at lower en- 
ergies. It is interpreted as a signature of a strong collective 
behavior of matter created in such collisions. The observa- 
tion of the development, with increasing multiplicity, of such 
slope in pp collisions is consistent with the hypothesis, that 
the larger the produced multiplicity, the more self-interacting 
and collective is the produced system. 

Thirdly, all the measured radii grow with event multiplic- 
ity, in each range separately. This is shown more clearly in 
Fig-El where we plot the radii as a function of {dNch/dr\)^^^^^ 
(For our pp data we use the (dA^^./i/drj)''''^^ |/Vr/,>i given in 
Tab. IB. R%i^ and Rf grow linearly with the cube root of 



charged particle multiplicity, for all kj ranges. Data, at both 
energies, follow the same scaling. For 7?,'f„, the situation is 
similar for medium ranges. The lowest ^x points show the 
strongest growth with multiplicity, while the highest hardly 
grows at all. That is the result of the strong change of the 
slope of kj dependence with multiplicity, noted in the discus- 
sion of Fig.nn 

Similar multiplicity scaling was observed in heavy-ion col- 
lisions at RHIC energies and below. In Fig.[T3]we compare 
our results to heavy-ion results from collision energies above 
15 AGeV. This is the first time that one can directly compare 
pp and heavy-ion radii at the same (dA^^./j/drj), as we mea- 
sure (d;Vc/j/dT|) comparable to the one in peripheral AuAu 
and CuCu collisions at RHIC. Since the value of the radius 
strongly depends on kj, we carefully selected the results to 
have the same average (^x) = 0.4 GeV/c. The picture at other 
^x's is qualitatively similar While both the heavy-ion and 
pp data scale linearly with (d;V(/,/dr|)'''^'^', the slope of the 
dependence is clearly different, for all directions. The pp re- 
sults are systematically below the heavy-ion ones at similar 
multiplicity; therefore, the "universal" multiplicity scaling jst] 
observed in heavy-ion collisions does not hold for pp colli- 
sions at = 0.9 and 7 TeV. The pp radii do scale linearly 
with multiplicity, but with a different slope. 

We speculate that the difference comes from a different way 
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that the two types of collisions arrive at similar multiplicity. 
To produce a large number of particles in pp collision one 
needs a particularly energetic elementary collision that pro- 
duces a lot of soft particles. The region where they are created 
is on the order of the incoming proton size and the growth of 
the size with multiplicity comes from further reinteraction be- 
tween particles after they are born. In contrast, in heavy-ion 
collision we have a combination of many elementary nucleon 
scatterings, each of them producing a relatively low multiplic- 
ity. But each scattering happens at a different space-time point 
and it is the distribution of these points that mostly determines 
the final observed size. In this picture, one would expect the 
heavy-ion sizes to be larger than the ones observed in at 
the same multiplicity. 



B. Systematic uncertainty 

The correlation function is, to the first order, independent 
of the single particle acceptance and efficiency. We performed 
the analysis independently for many samples of data, that nat- 
urally had single particle efficiencies different by up to 5%. 
We analyzed positive and negative pions separately, data at 
two magnetic field polarities, data from three different month- 
long "LHC periods", each of them having a slightly different 
detector setup. Two-particle correlations from all these ana- 
lyzes were consistent within statistical errors. 

We studied the effect of momentum resolution on the cor- 
relation peak with the MC simulation of our detector At this 
low px, below 1 GeV/c, the momentum resolution for tracks 
reconstructed in the TPC is below 1 %. This was confirmed by 
several methods, including the reconstruction of tracks from 
cosmic rays, and comparison of the reconstructed mass 
peak position with the expected value. The smearing of single 
particle momenta does result in the smearing of the correla- 
tion peak: it makes it appear smaller and wider. We estimated 
that this changes the reconstructed radius by 1 % for the fem- 
toscopic size of 1 fm; the effect grows to 4% for the size of 
2 fm, as it corresponds to a narrower correlation peak. 

In contrast to single particle acceptance, the femtoscopic 
correlation function is sensitive to the two-track reconstruc- 
tion effects, usually called "splitting" and "merging". The 
"splitting" occurs when one track is mistakenly reconstructed 
as two. Both tracks have then very close momenta. This re- 
sults in a sharp correlation peak at low relative momentum. 
We have seen such effects in the data and we took several 
steps to remove them. Firstly, the requirement that the track 
is simultaneously reconstructed in the TPC and ITS decreases 
splitting significantly. In addition, each cluster in the TPC is 
flagged as "shared" if it is used in the reconstruction of more 
than one track. The split tracks tend to produce pairs which 
share most of their clusters; therefore, we removed pairs that 
share more than 5% of the TPC clusters. We also look for 
configurations where a single track is split in two segments in 
the TPC, e.g. by the TPC central membrane or a TPC sec- 
tor boundary. Such segments should be correctly connected 
in the tracking procedure to form a single track if the detector 
calibration is perfect. However, in a few rare cases this does 



TABLE II. Systematic uncertainty coming from varying "mini-jet" 
background height/width by 10%/5% up/down. 



kj (GeV/c) Ai^C, % A^G^^ % Ai?G„„ % 



long 



(0.13,0.2) 
(0.2, 0.3) 
(0.3, 0.4) 
(0.4, 0.5) 
(0.5, 0.6) 
(0.6, 0.7) 



4 
4 
4 
7 
9 
13 



not happen and a split track can appear Such pairs would 
consist of two tracks that have a relatively small number of 
TPC clusters and they would rarely both have a cluster in the 
same TPC padrow. Therefore, we count, for each pair, the 
number of times that both tracks have a separate (non-shared) 
cluster in a TPC padrow. Pairs for which this number is low 
are removed. Both selections are applied in the same way to 
the signal and background distributions. As a consequence, 
the fake low-momentum pairs from splitting are almost com- 
pletely removed, and the remaining ones are concentrated in 
a very narrow relative momentum q range, corresponding es- 
sentially to the first correlation function bin. The inclusion of 
this bin has a negligible effect on the fitting result; hence, we 
do not assign any systematic error on the fitting values from 
these procedures. 

Another two track effect is merging, where two distinct 
tracks are reconstructed as one, due to finite detector space- 
point resolution. The ALICE detector was specifically de- 
signed to cope with the track densities expected in heavy- 
ion PbH-Pb collisions, which are expected to be orders of 
magnitude higher than the ones measured in pp collisions. 
More specifically, the ITS detector granularity as well as TPC 
tracking procedure, which allows for cluster sharing between 
tracks, make merging unlikely. We confirmed with detailed 
MC simulation of our detector setup that merging, if it appears 
at all, would only affect the correlation function in the lowest 
q bin, which means that it would not affect the measured radii. 

In summary, the systematic uncertainty on the raw measure- 
ment, the correlation functions itself, is small. 

The most significant systematic uncertainty on the extracted 
radii comes from the fact that we rely on the MC simulation 
of the "mini-jet" underlying event correlations. We fix the 
parameters of the B function in Eq. ^ by fitting it to the cor- 
relations obtained from the MC generated events. We con- 
firmed with the analysis of the non-identical 7t+7t^ pairs that 
our Monte-Carlos of choice, the Perugia-0 tune of the PYTHIA 
6 model, and the PHOJET model reproduce the height and the 
width of the "mini-jet" peak with an accuracy better than 10%, 
except the first multiplicity range where the differences go up 
to 20% for the highest fcx range. We performed the fits to the 
correlation function varying the parameters A/, and B/, of the 
B function by ±10%, and A„ by 5%. The fit values for the 
case when A/,, B/, are decreased and A,v is increased (corre- 
sponding to smaller "mini-jet" correlations) are systematically 
below the standard values. For larger "mini-jet" correlations 
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TABLE III. Systematic uncertainty coming from comparing the fit 
values with background obtained from PHOJET and PYTHIA simula- 
tions. 



(GeV/c) AR^„ 




long 


(0.13, 0.2) 


7 


4 


2 


(0.2, 0.3) 


1 


1 


4 


(0.3, 0.4) 


1 


1 


4 


(0.4, 0.5) 


7 


2 


4 


(0.5, 0.6) 


7 


3 


4 


(0.6, 0.7) 


10 


6 


7 



TABLE V. Systematic uncertainty coming from comparing the fits to 
two independent 3D correlation function representations. 



kj (GeV/c) AR 


■S« % 






(0.13,0.2) 


9 


5 


15 


(0.2, 0.3) 


9 


7 


7 


(0.3, 0.4) 


4 


2 


2 


(0.4, 0.5) 


6 


2 


4 


(0.5, 0.6) 


8 


3 


4 


(0.6, 0.7) 


18 


6 


12 



TABLE IV. Systematic uncertainty coming from varying the maxi- 
mum fit range. 



(GeV/c) ARg„ % ARf.,^ % ARg % 



(0.13, 0.2) 
(0.2, 0.3) 
(0.3, 0.4) 
(0.4, 0.5) 
(0.5, 0.6) 
(0.6, 0.7) 



3 
4 
7 
7 
7 
10 



they are systematically above. The resulting relative system- 
atic uncertainty on all radii is given in Table |ll] The error is 
independent of multiplicity, except for the first and last mul- 
tiplicity ranges, where it is higher by 50%. This error is fully 
correlated between multiplicity/A:T ranges. 

Independently, we performed the fits with the PHOJET gen- 
erator and fixed the parameters of B from them. The difference 
in the final fitted radii between PYTHIA and PHOJET back- 
ground is taken as another component of the systematic error, 
shown in Table Hill 

Another effect, visible in Fig.|9] is that the traditional Gaus- 
sian functional form does not describe the shape of the corre- 
lation perfectly. As a result, the extracted radius depends on 
the range used in fitting. Generally, the larger the fitting range, 
the smaller the radius. We fixed our maximum fitting range to 
1 .2 GeV, which is be sufficient to cover all correlation struc- 
tures seen in data. We estimate that the remaining systematic 
uncertainty coming from the fitting range is less than 5%. 

We always performed all fits separately to correlations for 
7t+7t+ and 71^71^ pairs. They are expected to give the same 
source size; therefore the difference between them is taken as 
an additional component of the systematic uncertainty. 

We used two independent representations of the 3D cor- 
relation functions: the "Cartesian" one uses standard 3- 
dimensional histograms to store the signal and the mixed 
background. The SH one uses sets of 1 -dimensional his- 
tograms to store the SH components plus one 3D histogram to 
store the covariances between them (see section llll Al for more 
details). The fitting of the two representations, even though 
it uses the same mathematical formula is different from 
the technical point of view. The SH procedure is more robust 
against holes in the acceptance [23], visible in our data, e.g. in 



Fig. [3] In an ideal case both procedures should produce iden- 
tical fit results; therefore, we take the difference between the 
radii obtained from the two procedures as an estimate of the 
systematic uncertainty incurred by the fitting procedure itself. 
The error is shown in Table IVl as a function of kj. The large 
error at low kj is coming from the fact that the two proce- 
dures are sensitive to the holes in the acceptance in a different 
way. It reflects the experimental fact that, in these kj ranges, 
pairs in certain kinematic regions are not measured; therefore, 
the femtoscopic radius cannot be obtained with better accu- 
racy. In the highest kj range the "mini-jet" underlying corre- 
lation is highest and broadest. If our simple phenomenological 
parametrization of it does not perfectly describe its behavior 
in full 3D space, it can affect differently a fit in Cartesian and 
SH representations. 

In summary, the combined systematic error is 10% for all 
and multiplicity ranges except the ones at the lower and 
upper edge. It is 20% for the lowest and highest kj and for the 
lowest and highest multiplicity range at each collision energy. 
It is also never smaller than 0. 1 fm. 



C. Non-Gaussian fits 

In the discussion of Fig.|9]we note that the measured cor- 
relation function is not perfectly reproduced by a 3D Gaus- 
sian fit. In our previous work II13I1 and in the work of the 
CMS collaboration lfl4ll it was noted that the shape of the 1- 
dimensional correlation in the Pair Rest Frame is better de- 
scribed by an exponential shape. Also, model studies ifisll 
suggest that pion production at these energies has large contri- 
bution from strongly decaying resonances. This is confirmed 
by the observation of significant resonance peaks in the 7i+7t^ 
correlation functions, seen e.g. in Fig. [8] Resonances de- 
cay after random time governed by the exponential decay law, 
which transforms into an exponential shape in space via the 
pair velocity. By definition pair velocity exists in the out and 
long direction, and vanishes in side. It is then reasonable to 
attempt to fit the correlation with a functional form other than 
a simple Gaussian, at least for the out and long components. 

If we keep the assumption that the emission function fac- 
torizes into the out, side and long directions, we can write a 
general form of the pair emission function: 



S{r)=So{rMr,)S,{n). 



(8) 
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We can independently change each component. We stress, 
however, that only for a Gaussian there is an analytically 
known correspondence between the pair separation distribu- 
tion S and single particle emission function Si. Two com- 
monly used forms of S are exponential and Lorentzian. They 
have the desired feature that the integration in Eq. (|2]i can be 
analytically carried out and produce a Lorentzian and expo- 
nential in C respectively. In order to select the proper com- 
bination of functional forms we seek guidance from models. 
They suggest that at least in the out and long direction the 
emission function is not Gaussian, and in some cases seems 
to be well described by a Lorentzian. We performed a study 
of all 27 combinations of the fitting functions for selected 
multiplicity/^T ranges. We found that universally the out cor- 
relation function was best described by an exponential, corre- 
sponding to Lorentzian emission function, which agrees with 
model expectations. In contrast, the side direction is equally 
well described by a Gaussian or a Lorentzian: we chose the 
former because the Lorentzian correlation function would cor- 
respond to exponential pair emission function with a sharp 
peak at 0. We deem this unlikely, given that the models do not 
produce such shapes. In long, the correlation function is not 
Gaussian; hence, we chose the exponential shape in C for the 
fit. In conclusion, we postulate that the source has the follow- 
ing shape: 



Sir) 



rexp 



7 

r: 



4R^ 



nE ^ 



(9) 



which corresponds to the following form of the femtoscopic 
part of the correlation function formula: 



Cf = l +A-exp 



R'^ a- 



'^side Hside 



(10) 

In Figs.fT4land[T5]we show an example of the exponential- 
Gaussian-exponential fit to the correlation functions at mul- 
tiplicity 23 < N^h < 29 and kj in (0.3,0.4) GeV/c. In the 
SH representation we see improvements over the Gaussian fit 
from Fig. |9] The behavior in CfJ at low q is now well de- 
scribed. In C, the "wiggle" in the correlation is also repro- 
duced - this is possible because the functional forms for the 
out and side directions are now different. In the Cartesian 
projections the improvement is also seen, however it is not 
illustrated as clearly as in the SH. 

We then proceed with the fitting of the full set of 72 cor- 
relation functions. The resulting fit parameters are summa- 
rized in Fig. [16] The quality of the fit (judged by the value of 
X^/Ndof) is better than for the 3D Gaussian fit. The A, param- 
eter is higher by up to 0.2, as compared to the pure Gaussian 
fit, reflecting the fact that the new functional form accounts for 
the pairs contributing to the narrow correlation peak at small 
q. The resulting exponential radii cannot be directly compared 
in magnitude to the Gaussian radii from other experiments. 
However all the features seen in dependencies of the Gaus- 
sian radii on multiplicity and kj are also visible here. This 
confirms that with a functional form that fits our correlation 
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FIG. 14. Exponential-Gaussian-exponential fit example for events 
with 23 < Nch < 29, pairs with 0.3 < Att < 0.4 GeV/c. SH represen- 
tation. Dotted line shows the femtoscopic component, dash-dotted 
shows the background, the dashed line shows the full fit. 
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FIG. 15. Exponential-Gaussian-exponential fit example for events 
with 23 < Nch < 29, pairs with 0.3 < ytx < 0.4 GeV/c. 1-dimensional 
projections of the Cartesian representation are shown, the other q 
components were integrated in the range — 0.16 GeV/c. 



function well (better than a 3D Gaussian) the physics message 
from the dependence of radii on multiplicity and kj remains 
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FIG. 16. Non-gaussian fit radii (see Eq. l llOt ) as a function of pair 
momentum kj for all multiplicity ranges and for two collision ener- 
gies. Panel a) shows ijf,,,, panel b) shows Rf/^fg, panel c) shows Rfg„^, 

panel d) shows Rout/^fide ratio. All points at given kj bin should be 
at the same value of kj, but we shifted them to improve visibility. 
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FIG. 17. 1 -dimensional Ri,„, radius for all multiplicity and kj ranges 
for the y/s = 0.9 TeV data. The points for different multiplicities 
were slightly shifted in kf for clarity. The systematic error, typically 
on the order of 10% is not shown I T 51. C losed and open stars show 
the previously published result from lOI for two ranges of the mul- 
tiplicity M. 



panel d) of Fig. [16] Again, the picture seen for the Gaussian 
radii is confirmed; the higher the multiplicity of the collision 
and the collision energy, the lower the value of the ratio. 



V. FITTING ID CORRELATIONS 

For completeness, we also repeated the 1 -dimensional 
study in Pair Rest Frame, using all the methods and fitting 
functions described in the previous work of ALICE ifisll . 
The 1 -dimensional correlation functions are fit with the stan- 
dard Gaussian form, modified with the approximate Bowler- 
Sinyukov formula to account for the Coulomb interaction be- 
tween charged pions: 



valid. The study of the fit functional form shows that the cor- 
relation does not have a Gaussian shape in out and long. 

The R'^jj^ from this fit should be equal to the Rfj^j^ from the 
3D Gaussian fit with two caveats. The first is the assumption 
that the emission function fully factorizes into separate func- 
tions for out, side, and long directions. In the fitting of the 
3D correlation functions the residual correlation between the 
value of the X parameter and the values of the radii is often 
observed. We noted already that the non-gaussian fit produces 
larger values of A,, so Rf^^^ could be affected. Nevertheless 
we observe very good agreement (within statistical errors for 
multiplicities above 16) between the values from both 
fits, giving us additional confidence that the underlying as- 
sumptions in our fit are valid. 

Similar conclusions can be drawn from the ratio of the 
Rfnit/R'^ide ^'^^ '■h^ more advanced functional form, shown in 



C{q,n.) = [{l-X)+XK{q,„,.) (l + exp{-Rl„.qf,j)] B{q„„), 

(11) 

where K is the Coulomb function averaged over a spherical 
source of the size 1 .0 fm, /?,„,, is the femtoscopic radius and 
B is the function describing the non-femtoscopic background. 
In Fig. [17] we plot the gaussian 1 -dimensional invariant ra- 
dius as a function of multiplicity and kj. The closed and open 
stars are the results from our earlier work, which are consis- 
tent with the more precise results from this analysis. The sys- 
tematic error is on the order of 10% and is now dominating 
the precision of the measurement. At ^ = 0.9 TeV we see 
that, for the lowest multiplicity, the radius is not falling with 
kj, while it develops a slope as one goes to higher multiplic- 
ity. The 1 -dimensional analysis is consistent with the 3D mea- 
surement — one needs to take into account that by going from 
LCMS (3D measurement) to PRF (ID measurement) because 
it is necessary to boost the out radius by pair velocity, which is 
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i/s = 7 TeV data. The points for different multiplicities were slightly 
shifted in kj for clarity. The systematic error, typically on the order 
of 10% is not shown Iil5i1 . 



defined by kj. Then, one averages the radii in three directions 
to obtain the ID /?,,„,. 

In Fig. [18] we show the same analysis performed for the 
^/s ~ 7 TeV data. The radii are again comparable at the 
same multiplicity/^:! range. In addition, as one goes to higher 
multiplicities, the kj dependence of /?,„,, is getting more pro- 
nounced. The results are again consistent with the 3D analy- 
sis. 



results in this work is that as one moves towards pp collisions 
producing high multiplicity of particles, similar collectivity 
develops. More experimental and theoretical information is 
needed to address this intriguing possibility. 

The upper range of multiplicities produced in p p collisions 
at y/s ~ 7 TeV is comparable to the multiplicities measured in 
peripheral heavy-ion collisions at RHIC. When plotted versus 
(d/Vc/i/dri)^'^^^ the radii in pp show linear scaling, but with 
different slope and offset than those observed in heavy-ion 
collisions. Therefore our observations violate the "universal" 
(dA?rf,/dri)''/^' scahng. This proves that the initial geometry 
of the collision does influence the final measured radii and 
must be taken into account in any scaling arguments. 

The analysis is complicated by the existence of the long- 
range underlying event correlations. We assume these are the 
"mini-jet" structures which are visible at values of pj as low 
as 0.5 GeV/c. The Monte-Carlo studies are consistent with 
such a hypothesis and are used to parametrize and take into 
account the influence of mini-jets on the fitted femtoscopic 
radii. Studies of the 7t+7t^ correlations are also consistent with 
such hypothesis. Nevertheless, the need to account for this 
effect remains the main source of the systematic error 

Finally, the detailed analysis of the correlation reveals that 
the three-dimensional Gaussian describes the measurement 
only approximately. A better shape, exponential-Gaussian- 
exponential, is postulated, based on Monte-Carlo studies, and 
is found to better agree with the data. The resulting radii and 
their behavior versus event multiplicity and pair momentum 
are fully consistent with the one obtained with the Gaussian 
approximation. 
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